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ABSTRACT 



We model the temperature distribution at the surface of a magnetized 
neutron star and study the effects on the observed X-ray spectra and light 
curves. General relativistic effects, i.e., red-shift and lensing, are fully taken 
into account. Atmospheric effects on the emitted spectral flux are not included: 
we only consider blackbody emission at the local effective temperature. In this 



first paper we restrict ourselves to dipolar fields. General features are studied 



and compared with the ROSAT data from the pulsars 0833-45 (Vela), 0656+14, 
0630+178 (Geminga), and 1055-52, the four cases for which there is strong 



evidence that thermal radiation from the stellar surface is detected. 



q . The composite spectra we obtain are not very different from a blackbody 

-M ■ spectrum at the star's effective temperature. We conclude that, as far 

3 i as blackbody spectra are considered, temperature estimates using single 

temperature models give results practically identical to our composite models. 
The change of the (composite blackbody) spectrum with the star's rotational 
phase is also not very large and may be unobservable in most cases. 

Gravitational lensing strongly suppresses the light curve pulsations. If a 
dipolar field is assumed, pulsed fractions comparable to the observed ones can 
only be obtained with stellar radii larger than what predicted by current models 
of neutron star structure, or with low stellar masses. Moreover, the shapes of the 
theoretical light curves with dipolar fields do not correspond to the observations. 
The use of magnetic spectra may rise the pulsed fraction sufficiently, but will 
certainly make the discrepancy with the light curve shapes worse: dipolar field 
are not sufficient to interpret the data. Many neutron star models with a meson 
condensate or hyperons predict very small radii, and hence very strong lensing, 
which will require highly non dipolar fields to be able to reproduce the observed 
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pulsed fractions, if possible at all: this may be a new tool to constrain the size 
of neutron stars. 

The pulsed fractions obtained in all our models increase with photon energy: 
the strong decrease observed in Geminga at energies 0.3 - 0.5 keV is definitely a 
genuine effect of the magnetic field on the spectrum in contradistinction to the 
magnetic effects on the surface temperature considered here. Thus, a detailed 
analysis of thermal emission from the four pulsars we consider will require both 
complex surface field configurations and the inclusion of magnetic effects in the 
atmosphere (i.e., on the emitted spectrum). 



Subject headings: dense matter stars: neutron stars: X-rays 



Submitted to The Astrophysical Journal 



-3- 

1. INTRODUCTION 

The detection of thermal emission from the surface of a neutron star is one of the 
'Holy Graal's of X-ray astronomy. At the first detection of X-rays from the direction of the 
Crab supernova remnant, surface thermal radiation from the neutron star already expected 
to be present in this remnant was proposed as the most probable source (Bowyer et al. 
1964), but the claim was soon disproved and even the most recent ROSAT observation of 
the Crab pulsar failed to detect any emission from the surface (Becker & Aschenbach 1993). 
A compilation of all Einstein observations (Seward & Wang 1988) listed ten neutron stars 
detected in soft X-rays, to which should now be added Geminga which at that time was not 
yet proved to be a neutron star but had been clearly seen by Einstein (Bignami, Caraveo 
& Lamb 1983). As of August 1993 there were thirteen confirmed detections of pulsars by 
ROSAT and six unconfirmed ones (Ogelman 1993). In most cases the radiation is probably 
of magnetospheric origin with some contamination from the surrounding synchrotron 
nebula, at least for the younger candidates, and surface thermal emission in some cases. 
With the sensitivity of the ROSAVs PSPC (Position Sensitive Proportional Counter) and 
long exposure times there is now strong spectral evidence that thermal radiation has been 
detected from four neutron stars (Ogelman 1993): PSR 0833-45 (Vela) (Ogelman, Finley 
& Zimmermann 1993), PSR 0656+14 (Finley, Ogelman & Kiziloglu 1992), PSR 0630+178 
(Geminga) (Halpern & Holt 1992), and PSR 1055-52 (Ogelman & Finley 1993). Moreover, 
these four objects show pulsed X-ray emission, three of them (except Geminga) are radio 
pulsars and three (except PSR 0656+14) have also been detected as 7-ray pulsars. 

Earlier Einstein observations had already provided some limited evidence for detection 
of surface thermal radiation in the cases of the Vela pulsar (Harnden et al. 1985), PSR 
0656+14 (Cordova et al. 1989), and Geminga (Halpern & Tytler 1988), while in the case 
of PSR 1055-52 (Cheng & Helfand 1983) thermal emission was considered as incompatible 
with the data. In later EXOSAT observations, with broad band spectroscopy only, thermal 
radiation was considered as the most reasonable origin of the detected soft X-rays in both 
cases of the Vela pulsar (Ogelman & Zimmermann 1989) and PSR 1055-52 (Brinkmann & 
Ogelman 1987). Because of its late discovery, as an X-ray source, in the Einstein data base, 
PSR 0656+14 has not been observed by EXOSAT. The EXOSAT observation of Geminga 
(Caraveo et al. 1984) did not give any new spectral information compared to the Einstein 
results. A review of the pre-ROSAT observational situation has been given by Ogelman 
(1991). 

The quality of the ROSAT data from these four nearby neutron stars presents a new 
chalenge for theorists to provide good models for their interpretation. The heretofore 
published analyses of these data have all assumed a unique surface temperature, with at 
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most a second thermal component coming from the small hot polar caps. We model here 
the temperature distribution at the surface of a magnetized neutron star and study its 
effects on the received spectra and light curves. The crustal magnetic field affects the heat 
transport in the layers beneath the surface and makes that regions of the star where the 
field is almost normal to the surface will be warmer than regions where the field is almost 
tangential to the surface. These temperature differences will give rise to modulation of 
the received flux at energies between 0.1 and 1 keV and are a natural explanation for 
the observed pulsations in the above mentioned four neutron stars. This long searched 
for detection of thermal radiation from the surface of neutron stars thus opens up a new 
window in the study of these objects. It has the potential to tell us about the structure 
of the surface magnetic field and give us new information about the size of these stars 
through gravitational lensing effects which obviously will be substantial. Our purpose in 
this first paper is to present the general physics involved and study the simplest case of a 
dipolar surface magnetic field. We have tried to present as clearly as possible the underlying 
physical ingredients as well as the method used, laying hopefully a clear groundstone for 
future improvements and/or modifications. 

General relativistic effects may be enormous in neutron stars and we take fully 
into account both gravitational red-shift and gravitational lensing. Magnetic fields are 
also affected by strong gravitational fields but we do not consider this effect: of critical 
importance for the surface temperature distribution is the angle between the surface's 
normal and the magnetic field, and this angle is practically unaffected by gravity (Ginzburg 
& Ozernoy 1964). The major effect of gravity on the surface magnetic field is to increase its 
strength at small radii: for a given dipolar field at infinity the surface strength is increased 
by 20-50% by gravity compared to the value it would have in flat space-time. The surface 
field strengths we consider should thus simply be somewhat reduced for comparison with 
values obtained for example from pulsar spin-down. Our results are however not very 
sensitive to changes of this size in the field strength. We do not include the effects of the 
magnetic field on the atmosphere where the emitted spectrum is determined. These effects 
are also substantial and will hopefully be included in future work. 

It has been proposed that the surface of neutron stars may be a magnetic solid 
(Ruderman 1974; Chen, Ruderman & Sutherland 1974). Even if improved calculations of 
the atomic lattice cohesive energy in very strong magnetic fields have indicated that this is 
not the case (Jones 1986), one should nevertheless keep in mind this possibility. It would 
obviously have dramatic effects on the emitted spectrum (Brinkmann 1980) and some, 
comparatively smaller, effect on the surface temperature (see Van Riper 1988 for simple 
estimates) Our model is based on magnetized envelope calculations which assume that 
the surface is not a magnetic solid: would this be uncorrect, our results should clearly be 
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reconsidered. 

We must also mention that the interpretation presented here assumes that the surface 
temperature is determined by the flow of heat from the star's interior and thus carries the 
imprint of the underlying crustal magnetic field. Another interpretation has been proposed 
(Halpern & Ruderman 1993) in which the hard X-rays emitted by the hot polar caps 
are scattered back onto the surface by the magnetospheric plasma: if this is the case the 
surface temperature has nothing to do with the properties of the magnetized envelope and 
our model is then obviously irrelevant. A composite model is also possible: the general 
temperature distribution may be determined by the heat flow from the interior and some 
regions heated by the back scattered hard X-rays. Moreover, absorption and scattering of 
radiation by the surrounding magnetospheric plasma may also be an important factor in 
reshaping the emitted flux (Halpern & Ruderman 1993). Which of these possibilities is 
actual can only be determined by studying each of them carefully. This paper is a first step 
in that direction. 

The structure of the paper is a follow. The ROSAT data are summarized in § |2] where 
we emphasize the features relevent to our purpose. In § |3] we describe the effects of the 
magnetic field on heat transport in the envelope and our model for the surface temperature 
distribution and in § f|] we present the method used to calculate the fluxes as observed by 
ROSAT. Our results are in § [5] and they are discussed in § |6], followed by our conclusions 
in§0 



2. THE ROSAT DATA 

We present in this section some characteristics of the ROSAT data for the four 
neutron stars we study along with some relevant informations obtained at other energies, 
summarized in Table [T]. Vela is the younger of the four and an accurate estimate of 
its surface temperature is of outmost importance since it could give us evidence for the 
occurence of fast neutrino emission in its early cooling history (Page & Applegate 1992). 
The other three are much older and their importance for cooling is distinct; they do provide 
us with insight about the state of their interior (Page 1994a). All four are X-ray pulsars, 
three of them (except Geminga) are also radio pulsars, three (except PSR 0656+14; see 
however Ma et al. 1993 and Brown Sz Hartmann 1993) are 7-ray pulsars and three (except 
PSR 1055-52) have an optical counterpart. Analysis of the X-ray emission from the Vela 
pulsar is delicate because of the strong emission from the surrounding synchrotron nebula 
which has to be substracted from the data to obtain the signal from the pulsar itself. The 
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other three X-ray sources are much cleaner than Vela: there is no evidence of extended 
emission in the ROSAT HRI observation of PSR 0656+14 (Anderson et al. 1993) and PSR 
1055-52 (Ogelman & Finley 1993). There is no HRI observation of Geminga to date but 
the PSPC data are compatible with a point source. 

All four objects show pulsed X-ray emission within the ROSAT energy range (0.08 
- 2.5 keV) and a two component spectrum, the soft component (roughly below 0.5 keV) 
being interpreted as surface thermal emission (Ogelman 1993). The corresponding pulsed 
fractions of these soft components are listed in Table |l[ Geminga's pulsed fraction below 
0.28 keV is larger than between 0.28 - 0.5 keV while the other three pulsars have a roughly 
constant pulsed fraction below 0.5 keV (a slight decrease seems to be present in PSR 
1055-52 below 0.5 keV). The phase of the peak in the light curve of these soft components 
is independent of the photon energy; the peaks are broad, the flux being above its mean 
value more than 50% of the time; the light curves are non sinusoidal. The hard component 
of Vela appears to be unpulsed: it may be real or just due to the fact that count rates 
are low in this part of the spectrum and the signal has been extracted from the nebula 
emission. In the case of PSR 0656+14 the count rate is too low to detect pulsation of the 
hard component but in the other two cases (Geminga and PSR 1055-52) pulsations are 
clearly seen and they are substantially off-phase compared to the soft component. The 
origin of the hard component of these spectra still remains ellusive due to the low count 
rates, fits with blackbody or power-law spectra giving equally good (or bad) results. It 
could be thermal emission from the hot polar caps or some magnetospheric process coming 
maybe, but not necessarily, from just above the polar caps. In this paper 'soft' and 'hard' 
X-rays will always refer to these two components within the ROSAT energy range. 

The effective temperatures T BB and interstellar column densities N^ B listed in Table [I] 
were obtained from single temperature blackbody fits. In the case of Vela, the quoted 
effective temperature and the observed flux imply a stellar radius of 3-4 km, i.e., only 
10% of the surface is emitting at this temperature: this may be due to very large surface 
temperature gradients or to a strong inadequacy of the blackbody spectrum. Use of more 
complex spectra from non magnetized He atmosphere (Romani 1987) or magnetized H 
atmosphere (Pavlov et al. 1994) give lower temperatures: 2.2 — 2.6 x 10 5 K (Finley et al. 
1992) and 6.6 - 7.4 x 10 5 K (Anderson et al. 1993) respectively instead of 8.6 - 9.4 x 10 5 
K with a blackbody spectrum in the case of PSR 0656+14. In the case of Geminga, 
magnetized hydrogen atmosphere spectra also give lower temperatures, 2 — 3 x 10 5 K 
(Meyer et al. 1994). 

The (dipolar) magnetic field strengths we cite are obtained the standard way (magnetic 
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dipole radiation braking) by 



B p = 2 x 10 12 (PP_ 15 ) 1/2 G 



where P is the period in seconds and P_is its derivative in units of 10~ 15 s s _1 . B p is thus 
the field strength at the magnetic pole on the star's surface. These values are of course only 
indicative, and general relativistic strengthening of the surface field adds to this uncertainty. 
Phenomenological analyses of the radio data allow to estimate the inclination angle a 
of the dipole with respect to the rotation axis and the angle ( between the observer's 
direction and the rotation axis, but the results are strongly model dependent in many cases 
(Miller & Hamilton 1993). PSR 1055-52 is generally considered as a typical example of 
radio inter-pulsar where emission from both magnetic poles is detected (a ~ ( ~ 90°); an 
alternative explanation for the interpulse is that the same magnetic pole is seen twice per 
rotation (a ~ ( ~ 0°). In the case of Geminga, modeling of the 7-ray emission within both 
the outer gap (Halpern & Ruderman 1993) and polar cap (Harding, Ozernoy & Usov 1993) 
models indicate that a ~ ( ~ 90°. The beaming of radio emission obviously requires a ~ ( 
for radio pulsars. 



3. STRUCTURE OF THE ENVELOPE 



For all neutron stars we are presently interested in, the interior is isothermal and 
a temperature gradient only appears in the upper layers of the outer crust (Nomoto & 
Tsuruta 1986, 1987). We will call these layers the envelope and fix its upper density at 
p = 10 10 gm/cm 3 . On top of the envelope resides the atmosphere: the total energy flux 
reaching it is determined within the underlying envelope, mostly by electron transport, 
and cannot be changed anymore, but its energy distribution, i.e., the emitted spectrum, is 
determined here. 



3.1. Envelope without magnetic field 

Detailed calculations of heat transport in non magnetized neutron star envelopes 
have been presented by Gudmundsson, Pethick & Epstein (1982, 1983) and Hernquist & 
Applegate (1984) These authors noticed that the study of the envelope can be separated 
from the study of the interior, since for neutron stars several months old the evolution time 



scale of the interior is much larger that the thermal time scale of the envelope, and the 
latter can be considered to be in a steady state. The structure of the envelope depends on 
the gravitational acceleration at the surface 

GM 

9s = , 2 

R\l-2GM/Rc 2 



and on the temperature at its base T&, i.e., at the envelope-interior boundary at density 
p = 10 10 gm cm -3 . Since the interior is isothermal its temperature is T int = TJ,. The envelope 
then determines the effective temperature T e related to the interior temperature T& quite 
accurately by 

/ rp4 \ 0.455 

T 68 = 1.3 MM (3) 

\9suJ 

(Gudmundsson, Pethick & Epstein 1982) where X n is X measured in 10™ cgs units. 
Since the interior is isothermal and, in absence of magnetic field, the heat transport 
in the envelope is isotropic the resulting surface temperature is uniform over the star: 

T s (e,<f)) = T e . 

The heat transport is due to electrons in the lower layers and to photons in the upper 
layers. The critical region, "sensitivity strip" , which determines the Tj, — T e relationship is 
the region just below the electron-photon transport transition where heat is thus transported 
mainly be electrons and the ions are in the liquid phase (Gudmundsson, Pethick & Epstein 
1982, 1983; Hernquist & Applegate 1984). When the surface temperature decreases below 
about ~ 3 • 10 5 K, the sensitivity strip is mostly in a region where the ions are partially 
ionized and there is no reliable calculation of the electron conductivity in this regime. 
Electron scattering by impurities also begins to contribute significantly; it depends on the 
impurity concentration, at most a poorly known quantity. Moreover the photon opacity 
has to be extrapolated significantly from calculated values. This makes that the surface 
temperatures cannot be reliably calculated and all such values (T s < 3 • 10 5 K) are only 
illustrative. The same remark apply a fortiori to the magnetic cases considered in the next 
subsection. 



3.2. Envelope with magnetic field 

In presence of a strong magnetic field the electron motion in directions perpendicular 
to the field becomes quantized, the electrons occupying Landau levels with spacings (non 
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relativistic case) 

hu B = 54^ = 11.3 k eV B 12 (4) 

m e c 

For typical fields present in neutron star envelopes electrons will occupy, at zero temperature, 
only the first Landau level up to densities of the order of 5 • 10 5 gm cm" 3 and then start to 
populate higher levels. Temperature can push electrons into higher levels only when 

k B T > %uj b . (5) 

The structure of the electron Fermi surface is thus strongly affected by the magnetic field 
almost across the whole sensitivity strip, and the transport properties will consequently 
depend on the field. Electron transport in directions perpendicular to the field is strongly 
suppressed while it is slightly enhanced parallelly to the field. Photon opacities are also 
affected but the Rosseland mean, the relevant quantity for heat transport, is only slightly 
anisotropic: anisotropy in heat transport thus comes mainly from electrons. A review of 
magnetized envelopes can be found in Yakovlev & Kaminker (1994). 

When compared to the non magnetic case, the surface temperature will be slightly 
higher if the magnetic field is radial (parallel transport) but much lower if the magnetic field 
is tangential to the surface (orthogonal transport). The envelope being at most a few tens 
of meters thick one can safely assume that the magnetic field is uniform across it. Detailed 
modelling of magnetized neutron star envelopes with uniform field have been presented by 
Hernquist (1985), Van Riper (1988) and Schaaf (1990a) for the case of parallel transport 
while the only published models for orthogonal transport are given by Schaaf (1990a). 
Hernquist's (1985) and Van Riper's (1988) results are in good agreement except at low 
temperature where Van Riper's inclusion of the (negative) ion Coulomb pressure (omitted 
by Hernquist) showed that the magnetic effects on the equation of state become dominant 
(Schaaf 's results are not in such good agreement due probably to the many approximations 
done). For the relevant field strengths, these magnetic T& — T s calculations can hence only 
be trusted at T s > 3 • 10 5 K, the same range of validity as for the non magnetic case. We 
will use Hernquist's (1985) results for parallel transport and Schaaf's (1990a) results for 
orthogonal transport. Figure [1] shows the resulting T& — T s relationships. 

In the particular cases of transport with a radial or a tangential field the problem 
is one dimensional since the heat flux can only be radial by symmetry. For an arbitrary 
orientation of the field with respect to the radial direction the problem of heat transport 
becomes two-dimensional. However, if the field orientation is sufficiently uniform on a 
length scale of a few tens of meters (= thickness of the envelope) the problem can locally 
be reduced to one dimension, at least in a first approximation. The surface temperature 
T s then depends only on the field strength B and the angle Q B between the field and the 
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radial direction (beside the dependence on the envelope's bottom temperature T b and the 
surface gravitational acceleration g s ). Greenstein & Hartke (1983) showed that one can 
then approximately write (keeping only the G# dependence explicit) 

T s (Q b )= X (Qb)xT s (Q b = 0) (6) 

where 

X (&b) = [cos 2 G B + xt sm 2 &b] ^ • (7) 

In Greenstein & Hartke's argument Xo — K±/K\\, the ratio of the thermal conductivities 
perpendicular and parallel to the magnetic field, which is assumed to be constant within 
the envelope. The Figure [l] gives the values of T s at G = 0° and 90° and we will use Equ. || 
and [7| to define T s at intermediate angles, thus taking 

_ T S (Q = tt/2) 
X ° = T s (Q = 0) ■ (8) 

Two main approximations are done in deriving Equ. |6| and |7|. The first one is that the 
radial temperature gradient dT/dr is much larger than the meridional temperature gradient 
dT/dx (x being a coordinate along the star's surface): T drops by about two orders of 
magnitude across the envelope whose thickness is only a few tens of meters, while it varies 
by less than a factor 10 (see Fig. [l]) along the surface on length scales of kilometers, thus 
dT/dr » dT/dx. The second approximation is that the thermal conductivities K\\ and K± 
are constant: Hernquist & Applegate (1984) showed that in fact, in the non magnetic case, 
the temperature profile T(p) within the outer part of the envelope where the conductivity 
is dominated by the photons (the region where most of the temperature drop occurs) 
does follow a trajectory of constant K. The validity of Equ. ^| and [7] can however only be 
checked by comparison with two-dimensional calculations. Schaaf (1990b) has performed 
such a calculation with relative success and many approximations which make it far from 
definitive, but his results are well reproduced by these simple formulas Equ. || and 0. 
Awaiting a more detailed study, Equ. § and |7| are the best that can be done. 



3.3. The atmosphere and the emerging spectrum 

The structure of the atmosphere is of critical importance for observations. The total 
flux permeating it is controlled by the underlying envelope but the spectral distribution of 
the flux (i.e., the emerging spectrum) is determined in this region. In the present work we 
will assume for simplicity that this emerging spectrum is a blackbody spectrum. Romani 
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(1987) took into account various possible chemical compositions of the atmosphere to 
calculate more realistic spectra in a zero magnetic field approximation. Due to the general 
uj~ 3 frequency dependence of the opacity k{uo) of a non magnetic fully ionized plasma, 
the specific flux is increased at high frequencies, reduced at low frequencies, and differs 
substantially from a blackbody spectrum. When partial ionization is taken into account the 
excess in the Wien part of the spectrum is reduced, the effect being most dramatic when 
metals are present. For most relevant temperatures, this part of the spectrum lies within 
the Einstein and ROSAT detector ranges. Inclusion of the magnetic effects is an enormous 
challenge and extensive calculations are in progress (Shibanov et al. 1992; Pavlov et al. 
1994). For a fully ionized magnetic plasma k b (uo) ~ uu~ l and the resulting spectrum is 
intermediate between the blackbody and the fully ionized non magnetic spectrum. Since 
the magnetic field increases substantially the electron binding energies, hydrogen absorption 
edges are present even at high temperature and alter the spectrum precisely within the 
detector ranges, pushing it toward the Planck spectrum or even below at low enough 
temperature (the details of these absorption edges are unfortunately blurred by the PSPC's 
low energy resolution). 

The blackbody approximation may finally not be as bad as it appeared originally in 
Romani's work. Moreover, the magnetic spectra depend on many parameters which have to 
be added to the parameter set of the models presented here and will make the analysis much 
more complicated. For these reasons, in this first work, we will consider only blackbody 
emission, characterized by the local effective temperature, and reserve more realistic spectra 
for future work. 



4. GENERATING LIGHT CURVES AND SPECTRA 



4.1. Flat Space-time 

Assume we have a neutron star at a distance D from the observer and we want to 
calculate the number of photons of energy E reaching a detector whose effective area for 
these photons is A(E). Choosing spherical coordinate (#, 0) on the star's surface with the 
coordinate axis pointing toward the detector, the number of photons emitted at the surface 
by an area d 2 A = R 2 sin Odd d<p, where R is the star's radius, is 

d xd u 

d 6 N = Phase Space Volume x Distribution Function = 2 — x n(E, T; . . .) (9) 

h 6 
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where the distribution function n may depend, beside the photon energy E and 
the temperature T of the emitting region d 2 A, on the local magnetic field, chemical 
composition, etc... For blackbody emission n = \j{e E l kBT — 1). Let's write 

d 3 x = cdt-d 2 Acos6 (10) 

(c is the speed of light and dt the time during which the received photons have been 
emitted) and 

7 o o , _ E 2 dE A(E) , x 

d 3 P = p 2 dp da p = — — ^ (ii) 



(dQ p = A(E)/D 2 is the detector's solid angle as seen from the star's surface and p = E/c 
the photon momentum), and integrate the emission over the portion of the star visible to the 
observer. Interstellar absorption at energy E is included through a factor exp(— Nh&{E)), 
where Nh is the effective hydrogen column density between the star and the observer and 
cr(E) the effective cross section (Morrison & McCammon 1983). We obtain thus 

d 2 N(E) = -^-^-< n (E; . . .) > e -^(E) E 2 A ^ df dE (12) 

c 2 h A D 2 



where 



<n(E;...)>= Clxdx F ' ^- n(E,T(6,<j>);...) (13) 

jo Jo 2n 



>o Jo 

is the distribution function averaged over the visible part of the star and x = sin6. 



4.2. Curved Space-time 

General relativistic effects are twofold here: red-shift and lensing. The red-shift makes 
that photons emitted at energy E are received at energy Eoo = Ee^ < E while photons 
emitted during a time dt are received in a time dtoo = e^^dt > dt where 



With this we can rewrite the photon flux in term of observed quantities as 

d 2 N(E 00 ) = -£- ^f < ra(e-* E^, . . .) > e -^( Eoc ) E 2 ^^ ^ ^ 

= M{E 00 )A{E 00 ) dtoo dE^ (15) 



-13- 

where 



Roo = Re -t> = R > R (is) 



1- 



2GM 



and the photon flux A/"(-Eoo) has dimension [photons cm -2 sec -1 keV -1 ]. 

The effect of gravitational lensing on photons emitted by hot spots on the surface of 
neutron stars has been studied in detail by Pechenick, Ftaclas & Cohen (1983). More than 
half of the star's surface is visible to an observer at infinity (Fig. |2|) and a photon emitted at 
collatitude 9 on the star's surface which reaches the observer must be emitted at an angle 6 
to the star's normal, defined through 

= e(x) = p- **> (17) 

m-)m) 2 -(l-2u)uV 



(1 



where 



x = sin 5. 



The impact parameter b of the photon trajectory is then given by b = x ■ Roq. In flat 
space-time 5 = 9 while here 5 < 9. The maximum collatitude visible to an observer 
at infinity corresponds to 5 = n/2, i.e., 9 max = 6(x = 1), and the star appears 
as a disc of radius Roo. For a given star mass, R^ is minimal at R = (3/2)Rg, 
with Rs = 2GM/c 2 = 2.95 km M/M & the star's Schwarzschild radius, and is 
R ao m = (3\/3/2) Rs — 2.QRs (thus, a 1.4 M neutron star radius cannot appear smaller 
than 10.74 km, for an actual value of 4.134 km). We show in Fig. ||] the value of the 
maximum lensing angle 9 max as a function of the star's radius R. 9 max reaches 180° at 
R/R s ^ 1.76, 360° at R/Rs — 1.5091, and goes to infinity when R/Rs tends to 1.5 
(r 7 = l.bRs is the so-called 'photon radius' at which a photon can be in a circular Keplerian 
orbit arround the star, if R < r~ 



i) 



To obtain the flux received by a detector the only change compared to the development 

of § |4.1.| is in Equ. [It] where cos# becomes cos S (the surface element d 2 A remains 

R 2 sin 9d6 dcj) by definition of the Schwarzschlid coordinates) and in the calculation of the 

detector's solid angle dfl p . The latter can be seen to be (see, e.g., Appendix I in Pechenick 

et al. 1983): 

'A(E)- 



x dx 
cos 5 sin 9 d6 



D 2 



(19) 



The term x dx in Equ. [13] is cos 9 sin 9 d9 which is exactly canceled by the term cos 8 sin 9 d9 
of Equ. ^ (since cos 9 — ► cos 5 in Equ. [TUp. The photon flux A/X-Eoo) is thus still given by 
Equ. [15] and O but now with x = sin 5. 
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The actual sizes of neutron stars are unknown and this lensing effect rises the hope of 
imposing some constraints on them. Theoretical predictions, for a 1.4 M star, range from 
6.5 km up to a maximum of 14 km. The smallest value (6.5 km) appears in models with 
either a kaon condensate (Thorsson, Prakash & Lattimer 1994) or hyperons (Pandharipande 
1971) and the largest values (14 km) are obtained in models with only nucleons and 
developed within the mean field theory. For models with only nucleons, the classical 
equation of state (EOS) of Friedman & Pandharipande (1981) and its improved versions by 
Wiringa, Fiks & Fabrocini (1988) give radii between 10.5 and 11.5 km at 1.4 M & . Very 
low mass neutron stars have much larger radii but, when restricting ourselves to masses 
arround 1.4 M & , 14 km should be considered as an upper limit in the models. 
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4.3. The detector's response 



The number of counts in the detector's output channel #i, Cts(i), is obtained from 
Af(E OQ ) with the detector's response matrix (r^): 

Cts{i) = At J2 r ij N( E j) SEj. (20) 

3 

For ROSAVs PSPC there are 256 output channels (i = 1, . . . , 256 = /, but channels 1 
- 7 have a zero effective area) and the incoming flux is split in 729 energy bands SEj 
(J = 1, . . . , 729 = J). At is the observing time. The r^ have dimension [cm -2 ] and the 
detector's effective area at energy Ej is A(Ej) = J2i r ij while the detector's effective area 
for channel i is A^ = J2j r ij- m an our calculations we use the 1992 March 19 release 
of the PSPC's response matrix. A perfect resolution detector would have I = J and 
Tij = Sij ■ A(Ei). 



4.4. Numerical method 



Given the star's mass M and radius R we define the magnetic dipole's strength B p 
and inclination angle a with respect to the rotation axis. The magnetic field is calculated 
at each point of a fixed grid on the star's surface (about one point per square degree) 
and, given the interior temperature T&, the surface temperature is calculated on this grid 
using Equ. ^| and |7J We then define a mobile grid on the surface, which covers the part of 
the star visible to the observer taking into account lensing, with usually one point per 3 2 
square degrees or one point per square degree for cases of very strong lensing. Given the 
observer's orientation ( with respect to the rotation axis, the mobile grid is rotated by steps 
of 6 degrees, at each step the temperature is calculated on this grid by interpolating from 
the values previously calculated on the fixed grid and the averaged distribution function 
(Equ. |13|) is calculated, usually for about one hundred values of the photon's energy E. 
Given the observer's distance D and the hydrogen column density Nh, the fluxes M{E) are 
calculated (Equ. [15]) and then passed through the detector's response matrix (Equ. |2(|) to 
give the counts which can be directly compared with the data. 
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5. RESULTS 

5.1. Surface Temperature Distributions 



With the method of § [3.2.| we can now calculate surface temperature distributions for 
given magnetic field configurations. Figure f| shows a typical example for a dipolar surface 
field, other dipole strengths and internal temperatures giving similar results. Figures |5|a) 
and b) give the percentage of the star's area with temperature inferior to a given value 
for various field intensities B p and internal temperatures T&. For a given TJ,, increasing B p 
results in a smaller proportion of the surface having a temperature close to the minimal 
temperature T min , while for a given B p the same effect is obtained by decreasing TJ,. Having 
the surface temperature distribution T s (9,4>) we can calculate the effective temperature 
T e =< T^(9,(j)) > 1 / 4 (< X > means surface average of X) and the average temperature 
T a =< T s {6, 0) > which will be used below. From Figures £| and | one sees that most of the 
surface is at a temperature close to the maximal temperature T max (reached arround the 
magnetic poles), but that the temperature drops to very low values on a small portion of 
the star. 



5.2. Spectra 

In this paper we use blackbody emission at each point on the surface, at the local 
temperature: our spectra are thus composite blackbody spectra. The effect of the non 
uniformity of the surface temperature on the spectrum turns out to be relatively small, the 
resulting spectra being very close to Planck spectra. Figure [5] shows a typical example, the 
same dipolar field as in Figure |J one sees that the blackbody spectrum at the effective 
temperature (dashed curves) is a reasonably good approximation to the 'exact' spectrum 
(continuous curves), particularly when the detector resolution is taken into account where 
both spectra are almost identical below 0.5 keV. Since a second spectral component is 
present above 0.5 keV, spectral fits to estimate the surface temperature are mostly sensitive 
to the soft band below 0.5 keV and consequently, as far as blackbody fits are concerned, 
single temperature fits are practically as good as our composite model to determine the 
effective temperature. Changing the field configuration and the observer's conditions 
(orientation, distance and column density) makes little difference. This result was already 
mentioned by Greenstein Sz Hartke (1983). The excess above 0.5 keV of the composite 
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spectrum compared to the blackbody spectrum at T e will though change the value of the 
parameters used to model the hard component. 

Phase dependent spectra are very close to the phase integrated spectrum, the difference 
being often below the detection limit with reasonable exposure times. Again, changing 
the dipole's strength and orientation, or the observer's conditions, makes little difference. 
If spectra from magnetic atmosphere models were used the results may be different since, 
when the temperature (and/or the surface field strength) varies, the shape of the spectrum 
changes while the blackbody spectrum's shape is invariant. Gravitational lensing has 
obviously no effect on the phase integrated spectrum, but will affect phase resolved spectra 
and make them even closer to the phase integrated spectrum. However light curves are 
more appropriate to study this effect. 



5.3. Light curves 

The properties of the observed light curves (see Table |TJ) that one should keep in mind 
here are the presence of one or two peaks, the phase of the peak(s) and the pulsed fraction, 
as well as their dependence on the photon energy. We define the pulsed fraction as: 

pf(-\ __ 2^ max ~ ^^ m »"/ (<-)i\ 

U iS mean 

where Cts max , Cts m i n and Cts mean are, respectively, the maximal, minimal, and mean 
count rates in a given detector channel (i=8,..., 256 for ROSATs PSPC) The observed 
pulsed fractions are of course delicate to measure and can only be defined for a range of 
channels, Pf{i\ — i?)- 

A surface dipolar field with two oppositely located warm regions on the star's surface 
gives light curves with one or two peaks, and singleness of the peak imply that the angle a 
between the magnetic dipole and the rotation axis and the angle ( between the observer's 
direction and the rotation axis must satisfy, in flat space-time, 

a + ( < 90° (22) 

to assure that only one magnetic pole is ever seen. Gravitational lensing will impose a 
more stringent condition but, on the other side, the statistical uncertainty in the observed 
counts makes that low amplitude peaks in the light-curve cannot be recognized. We show in 
Figure |7| seven examples of light curves with a = ( and the corresponding pulsed fractions. 
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We find, for single peak configurations, that the maximum pulsed fraction is obtained when 
a ~ £ ~ 45°, which is obvious from geometrical reasons. If a ~ ( > 45° the pulsed fraction 
is comparable or slightly lower than at 45° and a second peak appears. An important, 
simple, fact seen in all cases we have considered, from which only the most representative 
are presented in this paper, is that the phase of the peak(s) in the light curve is independent 
of the photon energy in the whole range 0.08 -2.5 keV, as can be seen for example in Fig. [7[ 

The lower the effective temperature T e the higher the pulsed fraction Pf which comes 
from the fact that at lower temperatures the difference between T max and T min is larger (see 
Fig. 0). A surprising result is that the highest magnetic fields do not give the largest pulsed 
fractions. As shown in Fig. |5|b, a larger portion of the surface is at low temperature at field 
strengths arround 3 ■ 10 11 G than at higher fields: it is thus at B p ~ 3 • 10 11 — 10 12 G that we 
obtain the highest pulsed fractions, and stronger fields have a pulsed fraction slightly lower. 

A general characteristic of pulsed fractions with dipolar fields and blackbody emission 
is that they rise with photon energy, slowly at low and high energy but with an intermediate 
region of faster increase whose exact location depends on the effective temperature but 
is always around 0.5 keV. The presence of this steeper rise around 0.5 keV is actually 
a product of the PSPC's response characteristics as can be seen from Fig. |S] where the 
pulsed fraction of the incoming flux (as it would be recorded by a detector with perfect 
energy resolution) is compared with PSPC-detected pulsed fractions. This feature will thus 
certainly be present even for more complex surface temperature distributions. 

If the detector has a very good energy resolution Pf is independent of interstellar 
absorption since its effect is just to reduce the incoming flux at energy E by a factor 
exp(— Nh o~{E)) which drops out when calculating Pf. However, when dealing with a 
detector like ROSATs PSPC, the finite energy resolution does affect the pulsed fraction 
by mixing in the same output channel photons with different energies and thus different 
absorptions. Since o~(E) ~ E~ 3 , the net effect of interstellar absorption is to suppress the 
soft photons compared to the harder ones and when increasing Nh the lowest PSPC's 
channels become proportionally more and more contaminated by hard photons. The result 
is a rise of the pulsed fraction in these lowest channels (Figure §). However, very large 
values of Nh are required to obtain a significant rise in Pf, and such high values are clearly 
excluded in the cases of the four pulsars we consider: they would dramatically reduce the 
flux at low energy, much below anything observed. 

The last point we consider is the effect of gravitational lensing. The effect on the 
spectrum is very small (at least for dipolar fields) and at the limit of observability due 
to the low count rates, but when adding the output of many channels into light curves 
their flattening can be clearly seen: the pulsed fraction is strongly dependent on the size 
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of the star, besides its much weaker dependence on the dipole's strength and the effective 
temperature. We show in Figure |9] a sequence of light curves where the star's radius is 
decreased from oo (i.e., no lensing effect) down to 6.25 km (i.e., R/Rs = 1-51), the star's 
mass being kept at 1.4 M Q . The corresponding pulsed fractions are shown in Fig. TO. The 
dipole and the observer are at 90° to the rotation axis and the dipole strength and interior 
temperature are choosen to maximalize the pulsations according to the above discussion. 
One immediately sees the dramatic effect of lensing which reduces Pf from about 25% 
(below 0.5 keV) in flat space-time to a half of this at a radius of 18 km (an unrealistically 
large value for a 1.4 M & neutron star) and to a tenth of this at a radius of 10 km. As 
another example, we took a series of 6 neutron stars with masses and radii of: (1.0 M Q , 
11.21 km), (1.2 M Q , 11.08 km), (1.4 M & , 10.89 km), (1.6 M , 10.55 km) and (1.8 M , 
9.76 km) (these radii correspond to the EOS of Friedman & Pandharipande 1981). With 
an internal temperature of 5 x 10 7 K, a dipole strength of 3 x 10 12 G, and Q = a = 90°, 
we obtained pulsed fractions, in the energy band 0.08-0.5 keV, of 12%, 8%, 4%, 1%, and 
1% respectively. There is also a tricky effect at radii between 8.85 km and 7.3 km: the 
peak of the light curve is 90° off-phase with respect to larger or smaller radii. At large 
radii (R > 8.85 km) the peaks correspond to one of the magnetic poles being seen face-on 
and when R decreases more of the cool region in the magnetic equator is visible at this 
phase, thus reducing the amplitude of the pulsations. When 7.3 < R < 8.85 km, the flux 
with the poles face-on is now lower than when they are seen sideway (phases 0.25 and 
0.75) where both warm polar regions are wholly seen because of the strong lensing; at still 
smaller radii, the second pole becomes visible in the face-on configurations and the peak is 
back at phases and 0.5. This latter effect of increase of the pulsed fraction at small radii 
when the whole surface of the star is visible and a hot region passes on the antipode with 
respect ot the observer, 'gravitational beaming', has already been described by Pechenick 
et al. (1983) and the off-phasing of the peak at slightly larger radii can also be seen in 
their figures. Orienting the magnetic dipole and the observer at 45° to the rotation axis, 
a = ( = 45°, gives very similar results: the light curves have the shapes shown in Fig. |7|, 
the pulsed fractions are observationally indinstinguishable from the ones of Fig. |TOJ (for the 
same temperature distribution) and the off-phasing of the peaks also happens at practically 
the same star's sizes. If we define a phase dependent effective temperature Tf from the 
average of T^(8, <ft) ■ cos 5 over the part of the star seen by the observer, i.e., an effective 
temperature that fits the phase dependent spectra, then, for the cases of Fig. |9|, Tf at the 
maximum exceeds Tf at the minimum by about 13% when there is no lensing, down to 6% 
at R = 18 km, 4.2% at R = 14 km and 0.9% at R = 10 km. 
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6. DISCUSSION 

The magnetic envelope calculations we use to model the surface temperature 
distribution tell us that the temperature of the warmer regions, where the field is radial 
(@b ~ 0), is several time higher than the temperature of the regions where the magnetic 
field is tangential (0# ~ 90°) as shown in Fig. [I]. There are thus very large temperature 
differences on the surface but the resulting effective temperature seen by an observer 
at infinity, who averages over the part of the star he sees, varies little during the star's 
rotation. Variations are of the order of 5-10% for observed pulsed fraction of 10-30% below 
0.5 keV and the observed spectrum varies thus little with the rotational phase. Magnetic 
effects on the emitted spectrum by themselves are important and it is possible that when 
included in our model they may produce more significant and interesting variations in the 
observed spectra. These will hopefully be included in future work and for the time being 
we will concentrate on the light curves. 

The first noticeable characteristic of the observed light curves for the four pulsars we 
study is their substantial pulsed fraction: the results of Fig. [10J show that these values, 
between 10 — 30% below 0.5 keV, could be obtained with simple dipolar fields. However the 
radii required, to avoid the flattening from gravitational lensing, are larger than the present 
theoretical predictions if the neutron star mass is assumed to be 1.4 M Q . Models with only 
nucleons in the core give the largest radii, and the most sophisticated models (Wiringa et 
al. 1988; Miither, Prakash & Ainsworth 1987) give radii between 10 - 12.5 km while other 
predictions barely reach 14 km (Glendenning 1985 e.g.). Models with meson condensates, 
hyperons or quarks give smaller radii (Thorsson et al. 1994, Glendenning 1985, Glendenning 
1992) For any of these values, the pulsed fractions in the energy range 0.08 - 0.5 keV, for 
whatever values of the parameters (B p , T&, a and (), are below 10%. Analyses of the Vela 
pulsar glitches do favor large radii (Link, Epstein & Van Riper 1992; Alpar et al. 1993) 
but values above 15 km (for a 1.4 M & neutron star) needed to give high enough Pf are 
not supported to date by any neutron star model. Older calculations of neutron stars with 
solid neutron cores did give radii of the order of 16 km (Pandharipande & Smith 1975), but 
the solidification of neutrons in neutron stars is not considered anymore as realistic. Taking 
neutron stars of lower masses raises the pulsed fraction, e.g., up to 12% for a 1.0 M star 
with radius of 11.2 km as seen in § |5.3.| using a realistic EOS, but still below 10% at masses 
above 1.2 M . So considering stars of mass lower than 1.4 M is not sufficient to obtain 
the high observed pulsed fractions unless a very stiff EOS is used. However, the values 
obtained are not much below the observed ones unless the star's radius is small. 

Spectra from magnetized hydrogen atmospheres give by themselves pulsed fractions of 
the order of 10% (Pavlov et al. 1994) for a dipolar field with a uniform surface temperature 
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and no lensing: when used with the appropriate temperature distributions it is reasonable 
to expect that they will increase Pf compared to the blackbody case, but the effect of 
the magnetic field considered in this paper will probably dominate (we obtain a 25% 
modulation without lensing). Using non magnetic hydrogen atmosphere spectra (Page, 
Zavlin, Pavlov & Shibanov, work in progress) one can get higher Pf, between 10-20% for 
standard masses and radii; the increase, compared to blackbody spectra, being due to the 
slight radial beaming produced by the limg darkening effect. 

Our model shows that, for a fixed magnetic field configuration, the difference between 
the maximum and minimum surface temperature and hence the observed pulsed fraction 
increases with decreasing temperature. It may be interesting to notice (see Table |XJ) that 
this trend is present within our very limited sample of four objects. It may just be a 
coincidence or may indicate the occurence of a general pattern of field evolution in which 
pulsars after a few thousands of years develop similar surface field configurations. 

In the cases of PSR 0833-45, 0656+14 and 1055-52 the observed pulsed fractions are 
weakly energy dependent below 0.5 keV, which is almost what we obtain. However, the 
slight rise in Pf in our results may be important since this does not correspond exactly 
to the observations: the observed pulsed fraction of PSR 1055-52 seems to be slightly 
decreasing with increasing energy (below 0.5 keV). In the case of Geminga, though, there 
is a significant decrease of Pf in the band 0.28 — 0.5 keV compared to its value below 
0.28 keV: this feature is absolutely impossible to reproduce in our models (we doubt that 
more complex field configurations will alter this conclusion). This may be evidence for the 
inadequacy of blackbody emission and the result of genuine atmospheric magnetic effects 
on the emitted spectrum. In case of a very high column density (see Fig. |[) we do obtain a 
decrease of Pf at low energy but it is much smaller than what is observed in Geminga and 
the column densities required, Ng > 10 21 cm -2 , are at least an order of magnitude too high. 

Even if we can somewhat reproduce the observed pulsed fractions we do not score 
well on the shape of the light curves in the three cases where they are single peaked (PSR 
0656+14, 0630+178 and 1055-52). A single-peaked light curve with a dipolar magnetic 
field restricts the geometry to a + ( < 90° and the resulting theoretical light curves have a 
long flat minimum and relatively narrow peaks (Fig. [7]): this is opposite to the observations 
where the peaks are broad and the light curves are above their mean more that 50% of the 
time. The use of magnetic spectra will change the shape of the light curves but probably in 
the 'wrong' direction: emission is beamed along the field and one can expect that the peaks 
will be narrower. We can safely conclude from this that the surface magnetic field of these 
three stars is certainly not dipolar. 

A simple way out of this problem is the addition of a small quadrupolar component 
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(Page 1994b) to broaden the peaks, keeping a + ( < 90°. That may be a reasonable solution 
for PSR 0656+14 but for PSR 1055-52 and 0630+178 the situation is probably different. 
Phenomenological description of the radio emission indicates (see Table |TJ) that a ~ ( ~ 30° 
or 8.2° for PSR 0656+14: if we accept this geometry then only one polar region is seen and 
the broadening of the soft X-ray light curve's peak should be achieved with an alteration 
of the dipolar field. However the small value of a and ( estimated by Lyne & Manchester 
(1988), 8.2°, yields pulsed fractions of at most 1% (Fig. |7|) with a dipolar field (no matter 
what the strength of the dipole or the temperature is); magnetic effects on the emitted 
spectrum will doubtly be sufficient to raise this to the observed 14%. The 30° — 35° 
inclination of Rankin (1993) or Malov (1990) give Pf ~ 8 — 10% and could work once 
magnetic spectra are used. A quadrupole component added to the dipole can also easily 
increase Pf though (Page 1994b). 

In the case of PSR 1055-52 and 0630+178 interpretation of the the radio and 7-ray 
emission, respectively, indicate that a ~ £ ~ 90°. A dipolar surface field will definitely 
give two peaks: if the dipolar component is orthogonal to the rotation axis there obviously 
must be strong higher order multipolar components, as already mentioned by Halpern & 
Ruderman (1993). This case will be considered in detail in a forthcoming paper (Page 
1994b). The alternative interpretation of PSR 1055-52 as being an almost aligned rotator 
(Malov 1990) was supported according to Malov (1989) by the absence of pulsation in the 
X-ray band in the Einstein observation (Cheng & Helfand, 1983), an argument invalidated 
by the ROSAT detection of pulsations. With a ~ ( ~ 30° (Malov 1990) we obtain a pulsed 
fraction just below 10% (Fig. 0) in the most optimistic case (B p = 10 12 G, low temperature 
and large radius of 16 km). Smaller radii will reduce Pf well below the observed 11% but, 
on the other hand, magnetic effect on the emitted sectrum will certainly be sufficient to 
give 11%. Thus the inclination angles proposed by Malov (1990) are probably able to give 
pulsed fractions as observed, but the asymmetric shape of the light curve again requires a 
non dipolar surface field. 

The light curves of PSR 0833-45 (Vela) are more delicate to interpret. The two ROSAT 
observations, separated by eight months, show two somewhat different light curves where 
several narrow subpeaks have moved in phase and amplitude. By considering the two 
observations together one can distinguish two peaks, the larger one about 210° wide and 
the smaller one about 100° wide with an amplitude about one third the amplitude of the 
large peak. However the dip which separates the small peak may be due to fluctuations 
and this small peak inexistent. An orthogonal configuration, a ~ ( ~ 90°, as suggested by 
analysis of the radio pulse, with a surface dipolar field, gives two soft X-ray peaks but the 
large difference in amplitude between the two observed peaks imply a distortion of the field 
from the pure dipole. An orientation with a ~ ( ~ 60° — 70° (as prefered, e.g., in outer 
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gap models of its 7-ray emission, Romani & Yadigaroglu 1994) would give better results, 
i.e., peaks of different amplitudes, but the peaks are definitely narrower than observed (and 
again inclusion of magnetic atmosphere effects will make the peaks still narrower). If the 
small peak is rejected then we are in the same situation as with PSR 0630+178 and 1055-52: 
there must be strong quadrupolar, or higher order, components in the magnetic field. 

Our discussion assumed that the contamination of the soft X-ray emission by the hard 
component is negligible: this is true when blackbody emission is used for both components, 
and also when magnetic atmosphere effects are included in the soft component (Meyer et 
al. 1994). Should the hard component contribute significantly to the soft band emission 
our conclusion may be affected. 

Several models of neutron stars with a pion or kaon condensate give very small radii 
(Maxwell & Weise 1976; Thorsson et al. 1994) and thus lead to very small modulations in 
the light curves, in sharp contradistinction to what is observed. We thus have a new tool 
to test these models, but one should not consider that our present results rule them out so 
far: surface magnetic field configurations more complex than the simple dipoles considered 
here, as well as magnetic effects on the emitted spectrum, may lead to stronger pulsations 
and the question needs more study (Page 1994b). If neutron stars are actually as small as 
these models with meson condensates claim, the off-phasing of the light curves with respect 
to the viewing (Fig. [J) at radii between 7.5 - 8.5 km must induce us to use extreme care in 
relating the soft X-ray peak phases with pulsations observed at other wavelengths. This 
range of radii includes the calculated values for a 1.4 M & neutron star in many models. 



7. CONCLUSIONS 

The simple model presented here, with dipolar fields, shows several features similar 
to the observed ones but with significant shortcomings. We obtain pulsed fractions which 
are lower that the observed ones, unless very large radii or very low masses are assumed. 
The shape of the light curves do not match the observations, indicating the presence of 
deviations from a purely dipolar field, the strength of the required correction depending on 
the assumed geometry and the size of the star. One can, however, reasonably expect that 
more complicated field structure will be able to reproduce the observed pulsed fractions and 
the shapes of the light curves: there is probably no need to invoke external heating of the 
surface and/or magnetospheric absorption (see Introduction) to reproduce the data, but 
this does not mean that these processes are not present. 'Energy dependent pulse shapes 
and phases are trying to tell us the whole story; we have to interpret them' (Ogelman 1993), 
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but deciphering the story will not be easy. A complete model will have to consider complex 
surface fields and also include the magnetic effects in the atmosphere. 

I am gratefull to J. H. Applegate, J. P. Halpern, Y. Shibanov and G. G. Pavlov for 
discussions, as well as A. Sarmiento, A. Serrano and E. Vazquez. My interest in this 
problem was triggered by an unpublished talk of S. Tsuruta (U. of Oklahoma, March 1991). 
This work was supported by a HEA-NASA grant NAGW 3075 at Columbia and by a 
DGAPA grant IN104092 at the UNAM. This is contribution number 350 of the IA-UNAM. 
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Table 1. Some observational properties 



PSR 0833-45 
(Vela) 



PSR 0656+14 



PSR 0630+178 
(Geminga) 



PSR 1055-52 



P (msec) 


89 




385 


t (yrs) 


1.2 x 10 4 


1 x 10 5 


T BB (K) b 


1.5- 1.6 x 10 6 


8.6- 9.4 x 10 5 


N BB (cm~ 2 ) b 


0.1 - 1.5 x 10 20 


0.8 - 1.2 x 10 20 


Soft component b 


Pf ~ 11% 
one or two peaks (?) 


Pf ~ 14% 
one peak 


Hard component 


Unpulsed 


? 


4>s - <t>H (deg) b 


n/a 


n/a 


Radio pulsar 


Yes (one pulse) 


Yes (one pulse) 


Optical counterpart 


Yes (pulsed) 


Yes (proposed) 


7-ray pulsar 


Yes (two pulses) 


No 


B P (G) 


6.8 x 10 12 


8.8 x 10 12 




( LM : 90 


( LM : 8.2 


a (deg) f 


I R: 90 


I R30 




[ M:35 


[ M35 


C (deg) f 


j LM : 83 

\ R:78 


< 


' LM : 8.2 
R: 30 

M:35 



Pf 



237 
3 x 10 5 

4.2 - 6.2 x 10 5 

1 - 2 x 10 20 

f 33% (0.08 - 0.28keV) 

[ 10% (0.28 - 0.5keV) 

one peak 

Pf ~ 35% 

~ 105 

No 
Yes (confirmed) 
Yes (two pulses) 



3.2 x 10 



90 



12 



90 



197 
5 x 10 5 

6.3 - 7.5 x 10 5 
2.0 - 4.5 x 10 20 

Pf ~ 11% 
one peak 



Pf ~ 20 



30% 



120 



Yes (two pulses) 

None to date 
Yes (one pulse) 

2.2 x 10 12 

LM : 74.7 

R90 

M30 

LM : 66.8 

R: 90 

M :30 



a The rows list respectively: the period P, spin down age r = P/2P, effective temperature from blackbody fit T BB and 

the associated hydrogen column density N^ B , pulsed fraction Pf and number of peaks of the soft X-ray component, pulsed 

fraction Pf of the hard X-ray component, phase difference between the peaks of the soft and hard X-ray components, 

detections at radio, optical and 7-ray wavelengths, the dipolar field strength (Equ^j) and the estimated angles between 

the magnetic dipole and the rotation axis (a) and between the observer's direction and the rotation axis (£). 

Values are from Ogelman et al. 1993 for Vela, Finley et al. 1992 for 0656+14, Halpern et al. 1993 for Geminga and 

Ogelman et al. 1993 for 1055-52. 
c Caraveo et al. 1994 

d Bignami et al. 1993 

e Fierro et al. 1993 

f Values are from Lyne et al. 1988 (LM), Rankin 1993 (R), Malov 1990 (M) and for Geminga: Halpern et al. 1993. 
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Fig. 1. — Tb — T s relationship for transport parallel and perpendicular to the field at various 
field strengths. Tb is the temperature at the bottom of the envelope whose density is fixed at 
10 10 gm cm~ 3 and T s is the 'surface' or effective temperature. The dependence on the star's 
mass and radius is entirely contained in the scaling factor g sl4: . The continuous line gives 
the non magnetic relationship. The high values of T s correspond to a radial magnetic field 
(parallel transport) and are taken from Hernquist (1985) while the low T s values correspond 
to a field tangential to the surface (orthogonal transport) and are from Schaaf (1990a). The 
values of T s lower than 10 55 K are linear extrapolations of the higher values and must be 
considered only as illustrative (see text). 

Fig. 2. — Geometry of gravitational lensing. 

Fig. 3. — Maximum lensing angle vs. neutron star's radius for a 1.4 M & star and vs. radius 
in unit of the Schwarzschild radius R s = 2GM/c 2 = 2.95 km M/M Q . 

Fig. 4. — Surface temperature distribution with a dipolar magnetic field with field strength 
B p — 3 ■ 10 12 G at the poles (marked by stars) and an internal temperature of 6 • 10 7 K 
(<7si4 = !)• The effective temperature is 6.27- 10 5 K and the average temperature is 5.8 • 10 5 
K. The representation is an area preserving mapping. 

Fig. 5. — Percentage of the star's area with temperature lower than a given value for dipolar 
fields. T m i n and T max are the minimum and maximum temperatures on the surface which 
can be read from Figure |l[ 

a) B p = 3 ■ 10 12 G with four different interior temperatures: 3 • 10 7 K (dash-dot), 6 ■ 10 7 K 
(dash-triple dot), 10 8 K (continuous) and 3 • 10 8 K (dash). 

b) T int = T b = 10 8 K with five different dipole strengths: 3 • 10 11 G (dash-dot), 10 12 G 
(dash-triple dot), 3 • 10 12 G (continuous), 6 • 10 12 G (dash) and 10 13 G (dot). 

These curves are independent of the mass and radius of the star. 

Fig. 6.— - Incoming (a) and detected (b) spectra for the temperature distribution shown in 
Figure ^. {D = 250 pc, Nh = 10 20 cm -2 and a = ( = 90°). The continuous lines show the 
total composite spectrum while the dotted lines are a blackbody spectrum at the average 
temperature and the dashed lines a blackbody spectrum at the effective temperature. The 
incoming spectra (a) are A/^-Eoo) and the detected spectra (b) are Cts(i)/At/6E, i.e., the 
incoming spectra passed through the response matrix. 
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Fig. 7.— Light curves (left) and pulsed fractions (right) for seven orientations of the 
observer and magnetic dipole: a = ( = 8°, 15°, 30°, 45°, 60°, 75° and 90°. In all cases the 
field is dipolar with strength B p = 10 12 G and the internal temperature T), = 5 x 10 7 K; 
T e = 5.2 x 10 5 K. 1.4 M neutron star with a radius of 16 km, N H = 10 20 cm~ 2 . The 
observer's distances are adjusted to avoid superposition of the light curves (the shape of the 
light curves and the value of the pulsed fraction are independent of the distance). 



Fig. 8.— Pulsed fraction for the configuration of Fig. [| and ^| but with different column 
densities N H as indicated. The continuous curve shows the pulsed fraction as would be 
obtained with a perfect detector, thus independent of the column density. 



Fig. 9.— Gravitational lensing of the light curves: orthogonal rotator. Sequence of light 
curves from a 1.4 M neutron star with decreasing radii from oo (i.e., no general relativistic 
effects) down to 6.25 km. Since the surface temperature, for a given field strength and internal 
temperature, is proportional to g sU (see Fig. P we have forced g s u = 1 in calculating the 
surface temperature distribution: all stars have exactly the same T s (9,(fi) and the effects 
shown are exclusively due to lensing and red-shift. We took T& = 5 x 10 7 K, B p = 10 12 G 
and a = ( = 90°. The distances of the stars are chosen such that the total received flux is 
the same for all stars and iV# = 10 20 cm -2 . Lowering of the intensity for decreasing radii 
is due to the red-shift which slowly drives the flux out of the PSPC's range. Notice that 
for radii between 8.85 km and 7.3 km the peaks are 90° off-phase with respect to the other 
cases. Phase corresponds to the dipole pointing toward the observer. 



Fig. 10.— Gravitational lensing of the pulsed fraction: same model as Fig. |9| with radii 
from oo down to 6.25 km. 



-31- 

Author's Address 

Dany Page: Instituto de Astronomia, U.N A.M., 

Apdo postal 70-264 

04510 MEXICO D.F., Mexico 

E-mail : PAGE@ASTROSCU.UNAM.MX 



Figure 2 




& 



Curved space-time trajectory 



Flat space-time trajectory 



GO 



Figure 1 



Figure 3 



6.5 



w 



o 



5.5 



1 1 


1 1 1 1 1 1 1 1 1 1 


' 1 | 1 1 1 . '. y 


_ _ . 


H = 10 14 G 


'is 


— . — 


H = 10 13 G 




— — 


■ H = 10 12 G ^ 
H = 10 11 G . jy 


^ y > 

■ ' / ; 




- h < io 10 <y^y 


.•'' / /'- 






■ - y yy 


• '// 




yyy 






/ y y - 


- 


.*x 


/ / - 






/ 


- ^v 




/ / '_ 


y 


y / y 

/ y . 


/ 


- 


/ y / 


- 


- 


/ / / 


. 


1 1 


i 1/ 1 i S i 1/ 1 i i 


■ i i i i i i 



5 

7 7.5 8 8.5 9 



360 



R (km), for M = 1 .4 M Q 

7.5 10 12.5 15 17.5 20 

t — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — n 



270 



X 

D 

E 



Q} 



180 - 



90 




1.5 2 2.5 3 3.5 4 4.5 5 



Log T. (K) 



R/Rc 



Figure 4 



7.4-10 K 



7-1 5 K 




6.27-1 (TK 



5.8-10 K 



3-1 5 K 



360 c 



180 £ 



Figure 5 a 



Figure 5b 



100 



3 10 12 G 



D 

< 50 




j i i i_ 



100 



Tint = 10° K 



D 

< 50 



mm 



max 




j i i i_ 



mm 



max 



Temperature 



Temperature 



I 
> 



o 

CO 
CO 



CM 
I 



E 
o 

CO 

c 
o 



o 

Q_ 



1 

0.1 
0.01 



10 
10 



-3 



-4 



10 



-5 



10 



-6 



2 10 



-7 



-8 



10 

10" 9 

-10 



10 



0.1 



Figure 6 



■ ■ I 


l 


r 


"^>\ a j 


r X 


^X 1 


r / 




-/ 


\\\ j 




A -J 


: 


\\\ j 


r 


* ' \ 1 


- 


\\\ j 


r 


' * \ a 


- 


\\\ j 


r 


' ' \ 1 


■ 


'AU 


r 


• * a 


■ ■ " 


■ '\1 



0.2 



0.5 




0.2 0.5 

Energy (keV) 



Figure 7a 



Figure 7b 



• l-H 

CO 

C 



CO 

c 




0.25 



0.2 - 



c 

o 0.15 

CD 

3 0.1 

PL, 



0.05 - 








Phase 



0.5 1 

Channel energy (keV) 



Figure 8 



a 
o 

^— » 
o 



7d 

CD 



0.25 



0.2 - 



0.15 - 



0.1 5 



0.05 - 







Perfect detector 




N^= 10 2 cm" 2 

N^ 10 2 cm" 2 

- N^= 10 2 cm" 2 

19 9 

N^ 10 cm" 2 



j L 



J L 



0.5 1 

Channel energy (keV) 



Figure 9 




0.5 1 1.5 

Phase 



CD 



0.5 



0.4 - 



O 0.3 

^— » 
o 



0.2 - 



0.1 - 







Figure 10 





0.5 1 

Channel energy (keV) 



